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Abstract 

Inclusive production of D*^ (2010) mesons in deep inelastic scattering has been 
measured with the ZEUS detector at HERA using an integrated luminosity of 
81.9 pb~^. The decay channel D*~^ — > D^7r~^ with D° K~7r~^ and correspond- 
ing antiparticle decay were used to identify D* mesons. Differential D* cross 
sections with 1.5 < < 1000 GeV^ and 0.02 < y < 0.7 in the kinematic region 
1.5 < Pt{D*) < 15 GeV and | 7]{D*) \ < 1.5 are compared to different QCD cal- 
culations incorporating different parameterisations of the parton densities in the 
proton. The data show sensitivity to the gluon distribution in the proton and 
are reasonably well described by next-to-leading-order QCD with the ZEUS NLO 
QCD fit used as the input parton density in the proton. The observed cross sec- 
tion is extrapolated to the full kinematic region in pt{D*) and ti{D*) in order 
to determine the open-charm contribution, F2'^{x,Q'^), to the proton structure 
function, F2. Since, at low Q^, the uncertainties of the data are comparable to 
those from the QCD fit, the measured differential cross sections in y and 
should be used in future fits to constrain the gluon density. 
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1 Introduction 



Charm quarks are produced copiously in deep inelastic scattering (DIS) at HERA. At 
sufficiently high photon virtuahties, , the production of charm quarks constitutes up to 
30% of the total cross section [TJ|2] . Previous measurements of D* cross sections [TJI21E1I11 
indicate that the production of charm quarks in DIS in the range \ < < 600 GeV^ 
is consistent with calculations in Quantum Chromodynamics (QCD) in which charm is 
produced through the boson-gluon-fusion (BGF) mechanism. This implies that the charm 
cross section is directly sensitive to the gluon density in the proton. 

In this paper, measurements of the D* cross section are presented with improved preci- 
sion and in a kinematic region extending to higher than the previous ZEUS results pp. 
Single differential cross sections have been measured as a function of and the Bjorken 
scaling variable, x. Cross sections have also been measured in two ranges as a func- 
tion of transverse momentum, pt{D*), and pseudorapidity, rj^D*), of the D* meson. The 
cross sections are compared to the predictions of leading-logarithmic Monte Carlo (MC) 
simulations and to a next-to-leading-order (NLO) QCD calculation using various parton 
density functions (PDFs) in the proton. In particular, the data are compared to calcu- 
lations using the recent ZEUS NLO QCD fit [5], in which the parton densities in the 
proton are parameterised by performing ffis to inclusive DIS measurements from ZEUS 
and fixed-target experiments. The cross-section measurements are used to extract the 
charm contribution, i^l"^, to the proton structure function, F2. 



2 Experimental set-up 

The analysis was performed with data taken from 1998 to 2000, when HERA collided 
electrons or positrons with energy = 27.5 GeV with protons of energy Ep = 920 GeV. 
The results are based on e~p and e~^p samples corresponding to integrated luminosities 
of 16.7 ± 0.3 pb~^ and 65.2 ± 1.5 pb~^, respectively.^ 

A detailed description of the ZEUS detector can be found elsewhere jH]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organised in nine superlayers covering the 



Hereafter, both electrons and positrons arc referred to as electrons, unless explicitly stated otherwise. 
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polar-angle^ region 15° < 6* < 164°. The transverse-momentum resolution for full-length 
tracks is (t{pt)/pt = 0.0058pT © 0.0065 © 0.0014/pt, with pt in GeV. 

The high-resolution uranium-scintillator calorimeter (GAL) consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a{E)/E = 0.18/ \/E for electrons 
and (t{E)/E = 0.35//E for hadrons, with E in GeV. 

Presamplers (PRES) [HI are mounted in front of FCAL, BCAL and RCAL. They consist 
of scintillator tiles which detect particles originating from showers in the material between 
the interaction point and the calorimeter. This information was used to correct the energy 
of the scattered electron. The position of electrons scattered close to the electron beam 
direction is determined by a scintillator strip detector (SRTD) jTUj. The SRTD signals 
resolve single minimum-ionising particles and provide a transverse position resolution of 
3 mm. 

The luminosity was measured from the rate of the bremsstrahlung process ep e'jp, 
where the photon was measured in a lead-scintillator calorimeter jTTj placed in the HERA 
tunnel at Z = — 107 m. 

A three-level trigger system was used to select events online [HIE]- At the third level, 
events with both a reconstructed D* candidate and a scattered-electron candidate were 
kept for further analysis. The efficiency of the online D* reconstruction, determined 
relative to an inclusive DIS trigger, was generally above 95%. 

3 Theoretical predictions 

A variety of models to describe charm production in DIS have been constructed, based 
on many theoretical ideas. A comparison of the data with these models is complicated 
by the need to produce predictions for the limited range of acceptance of the detector 
in pt{D*) and ri{D*). The calculation used in this paper to compare with the measured 
cross sections is based on NLO QCD as described in Section 13.11 Monte Carlo models 
also provide calculations in the measured kinematic region; those used are discussed in 
Section ing Predictions of other models are briefly discussed in Section ESI Most of these 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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models only predict the total cross sections and cannot therefore be directly compared 
with the current data. 

3.1 NLO QCD calculations 

The NLO predictions for cc cross sections were obtained using the HVQDIS program ^3] 
based on the so-called fixed-flavour-number scheme (FFNS). In this scheme, only light 
quarks {u, d, s) are included in the initial-state proton as partons whose distributions 
obey the DGLAP equations and the cc is produced via the BGF mechanism ^H] 
with NLO corrections ^H]- The presence of the two large scales, and m^, can spoil the 
convergence of the perturbative series because the neglected terms of orders higher than 
contain log(Q^/m^) factors which can become large. Therefore, the results of HVQDIS 
are expected to be most accurate at ~ and to become less reliable when ^ m^. 

The following inputs have been used to obtain the predictions for D* production at NLO 
using the program HVQDIS. The recent ZEUS NLO QCD global fit |2] to structure- 
function data was used as the parameterisation of the proton PDFs. This fit was re- 
peated P7] in the FFNS, in which the PDF has three active quark flavours in the proton, 
and Aq(^j~, is set to 0.363 GeV. In this fit, the mass of the charm quark was set to L35 GeV; 
the same mass was therefore used in the HVQDIS calculation of the predictions. The 
renormalisation and factorisation scales were set to /i = \/Q^ + '^ml for charm produc- 
tion both in the fit and in the HVQDIS calculation. The charm fragmentation to a D* 
is carried out using the Peterson function ^H]. The hadronisation fraction, /(c D*), 
taken from combined e"'"e~ measurements, was set to 0.235 ^H] and the Peterson param- 
eter, e, was set to 0.035 [201 • The production cross section for charmonium states at 
HERA is larger than in high-energy e"'"e~ collisions. The effect of J/i/j production on the 
hadronisation fraction was estimated from data PT|I^ to be about 2% and was neglected. 

As an alternative to the Peterson fragmentation function, corrections were applied to the 
partons in the NLO calculation using the Aroma MC program |2S| (see Section 13. 2p 
which uses the Lund string fragmentation [21], modified for heavy quarks according to 
Bowler and leading- logarithmic parton showers. This correction was applied on a 
bin-by-bin basis to the NLO calculation for each cross section measured, according to 
the formula rf(T(i:)*)NLO+MC = c?o-(cc)nlo ■ C'had where Chad = da{D*)Mc/da{cc)MC- The 
shapes of the differential cross sections calculated at the parton level of the Aroma model 
agreed reasonably well with those calculated from the HVQDIS program. The effect of 
the choice of hadronisation scheme is discussed in Sections IHl and ^| 

To estimate the contribution of beauty production, the NLO calculation and hadronisation 
from the MC were combined, using da{b — > -D*)nlo+mc = da{bb)^^i^o " C'had where Chad = 
da{b D*)Mc/da{bb)MC- The ZEUS NLO QCD fit was used as the proton PDF, so that 
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the mass used in this fit, = 4.3 GeV, was also used in the HVQDIS program and /x 
was set to -y/Q^ + 4m^. The hadronisation fraction, f{b D*), was set to 0.173 

An alternate way to describe charm production in QCD is the variable-flavour-number 
scheme (VFNS) |771 l28j. In these calculations, an attempt is made to treat the heavy 
quarks correctly for all Q^. Therefore, at low Q^, charm is produced dynamically through 
the BGF process as in the FFNS, whereas, at higher Q^, heavy-quark parton densities 
are introduced. The transition between the two extremes is treated in different ways by 
different authors j2Zll2H]- The ZEUS NLO QCD fit has been performed in this scheme 
using the formalism of Roberts and Thorne P^lOU] . Predictions from such calculations are, 
however, only available for the total charm cross section; no calculation of D* production 
in the measured kinematic range is available. 

3.2 Monte Carlo models of charm production 

The MC programs Aroma and Cascade [SJ were also compared with the measured 
differential cross sections. In the Aroma MC program, charm is produced via the BGF 
process. Higher-order QCD effects are simulated in the leading-logarithmic approxima- 
tion with initial- and final-state radiation obeying DGLAP evolution. The mass of the 
charm quark was set to 1.5 GeV and the proton PDF chosen was CTEQ5F3 [221 • The 
Cascade MC model takes a different approach to the generation of the hard sub-process, 
in which heavy-quark production is simulated in the framework of the semi-hard or kx- 
factorisation approach j221El]- The matrix element used in Cascade is the off-shell LO 
BGF process pHl35j . The Cascade initial-state radiation is based on CCFM evolu- 
tion jnn], which includes ln(l/a:) terms in the perturbative expansion in addition to the 
In terms used in DGLAP evolution. To simulate final-state radiation. Cascade uses 
Pythia 5.7 [H7]. The cross section is calculated by convoluting the off-shell BGF matrix 
element with the unintegrated gluon density of the proton obtained from the CCFM fit to 
the HERA F2 data |2H| with rric = 1.5 GeV. For both Aroma and Cascade, the Lund 
string model is used for the fragmentation into hadrons, and f{c^D*) was set to 0.235. 

3.3 Other predictions of charm production 

The extraction of F^"^ performed in this paper (see Section irUll is model dependent and 
comparisons of F2'^ to the predictions of models other than that used to produce it are 
not in general valid. Thus, only the FFNS model, which was used to extract Fg^'^, was 
compared to the data. 

Several models of charm production jSHl were compared in the x and range of the 
measurements in this paper. As most only predict total cross sections, the comparison 
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was performed for Fg^'^. All models show similar trends, with differences typically less than 
20%. Since the differences are smaller than the current precision of the D* cross-section 
measurements, these models are not considered further. 

4 Kinematic reconstruction and event selection 

The kinematic variables Q^, x and the fraction of the electron energy transferred to 
the proton in its rest frame, y, can be reconstructed using a variety of methods, whose 
accuracy depends on the variable of interest and its range: 

• for the electron method (specified with the subscript e), the measured energy and 
angle of the scattered lepton are used; 

• the double angle (DA) method jlO] relies on the angles of the scattered lepton and the 
hadronic energy flow; 

• the Jacquet-Blondel (JB) method [H] is based entirely on measurements of the hadronic 
system; 

• the E-method jl21 uses both the scattered-lepton energy and measurements of the 
hadronic system. 

The reconstruction of and x was performed using the S-method, since it has better 
resolution at low than the DA method. At high Q^, the E-method and the DA method 
are similar, and both have better resolution than the electron method. 

The events were selected jlSiP by the following cuts: 

• the scattered electron was identified using a neural- network procedure [H]. Its energy, 
-Eg', was required to be larger than 10 GeV; 

• Ve < 0.95; 

• 2/jB > 0.02; 

• 40 < (5 < 60 GeV, where 5 = Ei{l — cos6'j) and Ei is the energy of the calorimeter 
cell i. The sum runs over all cells; 

• a primary vertex position determined from the tracks fitted to the vertex in the range 
I ^vertex I < 50 cm; 

• the impact point (X, Y) of the scattered lepton on the RCAL must lie outside the 
region 26 x 14 cm^ centred on X = F = 0. 

The angle of the scattered lepton was determined using either its impact position on 
the CAL inner face or a reconstructed track in the CTD. The SRTD information was 
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used, when available. The energy of the scattered lepton was corrected using the PRES, 
with additional corrections for non-uniformity due to geometric effects caused by cell and 
module boundaries. The quantity S was calculated from a combination of CAL clusters 
and tracks measured in the CTD. The contribution to 6 from the scattered lepton was 
evaluated separately after all corrections were applied as described above. 

The selected kinematic region was 1.5 < Q"^ < 1000 GeV^ and 0.02 < y < 0.7. 

5 Selection of D* candidates 

The D* mesons were identified using the decay channel D*~^ D^^^t with the subsequent 
decay K~ti^ and the corresponding antiparticle decay, where 7r+ refers to a low- 

momentum ("slow") pion accompanying the D^. 

Charged tracks measured by the CTD and assigned to the primary event vertex were 
selected. The transverse momentum was required to be greater than 0.12 GeV. Each 
track was required to reach at least the third superlayer of the CTD. These restrictions 
ensured that the track acceptance and momentum resolution were high. Tracks in the 
CTD with opposite charges and transverse momenta pt > 0.4 GeV were combined in pairs 
to form candidates. The tracks were alternately assigned the masses of a kaon and 
a pion and the invariant mass of the pair, Mk-k, was found. Each additional track, with 
charge opposite to that of the kaon track, was assigned the pion mass and combined with 
the Z)°-meson candidate to form a D* candidate. 

The signal regions for the reconstructed masses, M{D^) and AM = {Mk-k-ks ~ Mktt), 
were 1.80 < M{D°) < 1.92 GeV and 0.143 < AM < 0.148 GeV, respectively. To allow 
the background to be determined, candidates with wrong-sign combinations, in which 
both tracks forming the candidates have the same charge and the third track has the 
opposite charge, were also retained. The same kinematic restrictions were applied as for 
those D° candidates with correct-charge combinations. 

The kinematic region for D* candidates was 1.5 < pt{D*) < 15 GeV and \ri{D*)\ < 1.5. 
Figure n shows the AM distribution for the D* candidates together with the background 
from the wrong-charge combinations. The fit to the distribution has the form 

F = pi- exp (^-0.5 ■ x^'^iT^ j + Pi ■ (AM - m^)^^ 

where x = |(AM — P2)/P3\, Pi — Ps are free parameters and is the pion mass. The 
"modified" Gaussian was used to fit the mass peak since it gave a better value than 
the conventional Gaussian form for a MC sample of D* mesons. The fit gives a peak 
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at 145.49 ± 0.02(stat.) MeV compared with the PDG value of 145.421 ± 0.010 MeV gg. 
The measured peak position differs from the PDG value. However, it was not corrected 
for detector effects and the systematic uncertainty was not determined. The fitted width 
of 0.61 ± 0.02 MeV is consistent with the experimental resolution. Consistent results were 
also found for the e~^p and e~p data separately. For the range 0.143 < AM < 0.148 GeV, 
a clear signal of D'^ candidates is also shown in Figure ^ 

The number of D* candidates determined in the two signal regions and after subtracting 
the background estimated from the wrong-charge sample was 5545 ± 129. The normalisa- 
tion factor of the wrong-charge sample was determined as the ratio of events with correct- 
charge combinations to wrong-charge combinations in the region 150 < AM < 165 MeV. 
This factor is compatible with unity for both e~p and e~^p data. The normalisation factors 
were determined for each bin in order to calculate the differential cross sections using the 
background-subtraction method. 

6 Acceptance corrections 

The acceptances were calculated using the Rapgap 2.08 jlHI and Herwig 6.1 |17] MC 
models. The Rapgap MC model was interfaced with Heracles 4.6.1 HH] in order to 
incorporate first-order electroweak corrections. The generated events were then passed 
through a full simulation of the detector, separately for e~p and e~^p running, using Geant 
3.13 |in] and processed and selected with the same programs as used for the data. 

The MC models were used to produce charm by the BGF process only. The GRV94- 
LO jnOl PDF for the proton was used, and the charm-quark mass was set to 1.5 GeV. 
The Herwig MC contains leading-logarithmic parton showers whereas for Rapgap MC, 
the colour-dipole model jSI] as implemented in Ariadne 4.03 was used to simu- 
late QCD radiation. Charm fragmentation is implemented using either the Lund string 
fragmentation (Rapgap) or a cluster fragmentation [22] model (Herwig). 

Figure 121 shows distributions of DIS variables for D* events (after background subtraction) 
for data compared to detector-level Rapgap predictions. The distributions, which are 
normalised to unit shown separately for two intervals: 1.5 < < 1000 GeV^ 

and 40 < < 1000 GeV^. The Rapgap predictions are in good agreement with the 
data distributions for both the scattered-lepton and hadronic variables. The description 
is similarly good for the two ranges. This good description gives confidence in the 
use of the Rapgap MC to correct the data for detector effects. The Herwig MC gives 
a similarly good representation of the data (not shown) and is used to estimate the 
systematic uncertainty, arising from the model in the correction procedure, as described 
in Section ISl 
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The cross sections for a given observable Y were determined using 

da_ _ N 

dY ~ A-C- B-AY' 

where is the number of D* events in a bin of size Ay, A is the acceptance (which takes 
into account migrations, efficiencies and QED radiative effects for that bin) and C is the 
integrated luminosity. The product, B, of the appropriate branching ratios for the D* 
and L>° was set to (2.57 ± 0.06)% jlSj. 

7 D* rates in e~p and e^p interactions 

The D* production rate, r = N/C, in the e~p data set is systematically higher than 
that in the e~^p data set. This difference increases with Q^; for example, the ratio 
of the rates, r^'P/r^^p^ is equal to 1.12 ± 0.06 for 1.5 < < 1000 GeV^, while for 
40 < < 1000 GeV^ it is 1.67±0.21 (only statistical errors are given). Such a difference 
in production cross sections is not expected from known physics processes. 

A detailed study was performed to understand whether any instrumental effects could 
account for the difference between the two data sets. No such effect was seen in inclusive 
DIS where the ratio of e~p to e^p rates is consistent with unity. The rate for the wrong- 
charge background under the D* mass peak in e~p data agreed well with the wrong- 
charge rate in e^p data. For example, for > 40 GeV^, where the largest difference 
exists, the ratio of the rates for wrong-charge track combinations in e~p and e~^p data is 
0.95 ±0.09. For both e~p and e'^p interactions, the number of D*^ mesons was consistent 
with the number of D*~ for the entire range studied. Different reconstruction methods, 
cuts, background-subtraction methods and the time dependence of the difference were 
also investigated. None of these checks gave an indication of the source of the observed 
difference between the D* rates in e~p and e~^p for > 40 GeV^. The cross sections were 
measured separately for e~p and e~^p data and are discussed in Sectional The difference 
in observed rate is assumed to be a statistical fluctuation and the two sets of data were 
combined for the flnal results. 

8 Experimental and theoretical uncertainties 
8.1 Experimental uncertainties 

The systematic uncertainties of the measured cross sections were determined by changing 
the selection cuts or the analysis procedure in turn and repeating the extraction of the 
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cross sections jSS|- The following systematic studies have been carried out (the resulting 
uncertainty on the total cross section is given in parentheses): 

• event reconstruction and selection The following systematic checks were per- 
formed for this category: 

— the cut on i/e was changed to i/e < 0.90; 

— the cut on yj^ was changed to i/jb > 0.03; 

— the cut on 6 was changed to 42 < 5 < 57 GeV; 

— the cut on the l^'vertcxl was changed to | ^vertex I < 45 cm; 

— the cut on Ef,' was changed to Ef.' > 11 GeV; 

— the cut on the position of the scattered lepton in the RCAL was increased by 1 cm; 

— the electron method was used, except for cases when the scattered-lepton track 
was reconstructed by the CTD. In the latter case, the DA method, which has the 
best resolution at high Q^, was used; 

— the energy of the scattered electron was raised and lowered by 1% in the MC only, 
to account for the uncertainty in the GAL energy scale; 

— the energy of the hadronic system was raised and lowered by 3% in the MG only, 
to account for the uncertainty in the hadronic GAL energy scale; 

— the reconstructed SRTD hit position was shifted by ±2 mm to account for the 
uncertainty in the SRTD-RGAL alignment. 

• uncertainties related to the D* reconstruction (li'e%). The following systematic checks 
were performed for this category: 

— tracks were required to have | ^7 |< 1.75, in addition to the requirement on the 
number of superlayers; 

— the cut on the minimum transverse momentum for the tt and K candidates was 
raised and lowered by 0.1 GeV; 

— the cut on the minimum transverse momentum for the tTs was raised and lowered 
by 0.02 GeV; 

— the signal region for the M(D°) was widened and narrowed symmetrically around 
the centre by 0.01 GeV; 

— the signal region for the AM was widened symmetrically around the centre by 
0.003 GeV. 

• the acceptance was determined using Herwig instead of Rapgap (—2.7%); 

• the uncertainty in the luminosity measurement (2.2%). 
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The cross section obtained using the fit was in good agreement with that obtained by 
subtracting the background using the wrong-charge candidates. These estimations were 
also made in each bin in which the differential cross sections were measured. The overall 
systematic uncertainty was determined by adding the above uncertainties in quadrature. 
The normalisation uncertainties due to the luminosity-measurement error, and those due 
to the D* and decay branching ratios of 2.5% jlH], were not included in the systematic 
uncertainties for the differential cross sections. 



8.2 Theoretical uncertainties 

The NLO QCD predictions for D* production are affected by the systematic uncertainties 
listed below. Typical values for the systematic uncertainty are quoted for the total cross 
section: 

• the proton PDF. The CTEQ5F3 and GRV98-HO [HH PDFs were used to check the 
sensitivity of the predictions to different parameterisations of the gluon density in the 
proton. The appropriate masses used in the fit to determine the PDF were also used 
in HVQDIS, i.e. 1.3 GeV for CTEQ5F3 and 1.4 GeV for GRV98-HO. The change in 
the cross section was +2.0% using CTEQ5F3 and -16% using GRV98-HO; 

• the mass of the charm quark (191%) . The charm mass was changed consistently in 
the PDF fit and in HVQDIS by +0.15 GeV. The largest effect was at low Pt{,D*); 

• the renormalisation and factorisation scale, jj, (li%) . The scale was changed by a factor 
of 0.5 and 2; another scale, 2mc, was also used The maximum of \J + 
and 2mc as a function of was taken as the scale to estimate the upward uncertainty; 

• the ZEUS PDF uncertainties propagated from the experimental uncertainties of the 
fitted data (±5%). The change in the cross section was independent of the kinematic 
region; 

• uncertainty in the fragmentation (l4%) . The parameter e in the Peterson fragmenta- 
tion function was changed by ±0.015. 

The first source of systematic uncertainty is shown separately in the figures. The last 
four were added in quadrature and displayed as a band in the figures. An additional 
normalisation uncertainty of 3% |H| on the hadronisation fraction /(c — > D*) is not 
shown. 
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9 Cross-section measurements 



9.1 Visible cross sections 

The overall acceptance after applying the selection criteria described in Sections |4] and El 
for 1.5 < g2 < 1000 GeV^, 0.02 < y < 0.7, 1.5 < Pt{D*) < 15 GeV and |r/(D*)| < 1.5 
calculated with Rapgap is 31%, both for e~p and e^p data. The total cross sections in 
the same region are: 

a{e-p e-D*X) = 9.37 ± 0.44(stat.)l{]:|(syst.) ± 0.23(BR) nb; 

(T(e+p ^ e^D*X) = 8.20 ± 0.22(stat.)+°:^^(syst.) ± 0.20(BR) nb, 

where the final uncertainty arises from the uncertainty on the branching ratios for the D* 
and D^. The D* cross section for e'^p data is consistent with the previously published 
result P obtained at a proton beam energy of 820 GeV. According to HVQDIS, a 5% 
increase in the D* cross section is expected when the proton energy increases from 820 to 
920 GeV. 

The cross section obtained from the combined sample is: 

a{e^p e^D*X) = 8.44 ± 0.20(stat.);!]J^(syst.) ± 0.21(BR) nb. 

The prediction from the HVQDIS program is 8.41 'tl'gl nb, in good agreement with the 
data. The uncertainty in the HVQDIS prediction arises from the sources discussed in 
Section IH21 (excluding that from using a different proton PDF) and is about 2.5 times the 
size of the uncertainty in the measurement. A contribution to the total cross sections arises 
from D* mesons produced in bb events. The D* cross section arising from bb production 
was estimated, as described in Sectional to be 0.17 nb for > 1.5 GeV^. The measured 
differential cross sections include a component from beauty production. Therefore, all 
NLO predictions include a bb contribution calculated in each bin. For the extraction of 
Fl'^, the predicted value of bb production was subtracted from the data. 

9.2 Differential cross-section measurements 

The differential D* cross sections as a function of Q^, x, pt{D*) and t]{D*) for the com- 
bined e~p and e^p data samples are shown in Fig. 01 and given in Table Q The cross 
sections in and x both fall by about four orders of magnitude in the measured region. 
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The cross-section da / dpTiD*) falls by two orders of magnitude with increasing pxiD*). 
The cross-section da/dri{D*) rises with increasing ri{D*). The ratio of the e~p and e~^p 
cross sections, also shown in Fig. El and given in Tabled tends to increase with increasing 
and x. Neither the NLO calculations nor the MCs based on LO matrix elements and 
parton showers depend on the charge of the lepton in ep interactions. 

The data in Fig. El are compared with predictions from the MC generators Aroma and 
Cascade. The prediction from Aroma is generally below the data, particularly at low 
and medium to high pt{D*). In contrast, the prediction from Cascade agrees at 
low Q^, but generally lies above the data. Both MC predictions describe the shapes of 
the cross-sections da/dx and da/dr]{D*) reasonably well. The uncertainties in these MC 
predictions are difficult to estimate and may be large. 

In Fig. m the same data are compared with the NLO calculation implemented in the 
HVQDIS program. The predictions used the default parameter settings as discussed in 
Section 121 with the uncertainties described in Section IH^ Predictions using an alternate 
PDF, CTEQ5F3, and an alternate hadronisation scheme, from Aroma, are also shown. 
The differences between the predictions, which are comparable to the uncertainties in 
the data, demonstrate the sensitivity of this measurement to the gluon distribution in the 
proton. The ratio of data to theory is displayed for each variable. For the cross sections as 
a function of and x, the NLO predictions give a reasonable description of the data over 
four orders of magnitude in the cross section. For da/dQ"^, the description of the data is 
similar over the whole range in Q^, even though HVQDIS is expected to be most accurate 
when ~ m^. The NLO calculation does, however, exhibit a somewhat different shape, 
particularly for da/dx, where the NLO is below the data at low x and above the data at 
high X. The predictions using CTEQ5F3 instead of the ZEUS NLO fit, or using Aroma 
for the hadronisation instead of the Peterson function, give better agreement with the 
data for the cross-section da/dx. 

The cross sections as a function of pt{D*) and ri{D*) are also reasonably well described 
by the NLO calculation. The prediction using the ZEUS NLO QCD fit gives a bet- 
ter description than that using CTEQ5F3 (and also better than the prediction using 
GRV98-HO, not shown), especially for the cross-section da/dr]{D*). A better description 
of da/dr]{D*) is also achieved [HS] by using Aroma for the hadronisation, although, in 
this case, da / dpT{D*) is not so well described. It should be noted that previous publica- 
tions P1I2] revealed discrepancies in the forward ri{D*) direction. This region can now be 
reasonably well described by a recent fit to the proton PDF as shown in Fig. Efd). The 
data presented here are practically independent of the data used in the ZEUS NLO PDF 
fit to inclusive DIS data. Further refinement of NLO QCD fits and even the use of these 
data in future fits may achieve a better description. 

Cross sections as a function of ri{D*) and pt{D*) were also measured for > 40 GeV^. 
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The combined e~p and e~^p data samples are given in Table|21and shown in Fig. compared 
with the HVQDIS predictions. Although the HVQDIS calculation is not thought to be 
applicable at high Q^, the data are well described. The high-Q^ region is also where the 
difference in e~p and e~^p data is most pronounced; the ratios of the cross sections are 
given in Table |21 



10 Extraction of Fg^^ 

The open-charm contribution, ^2^^^, to the proton structure-function F2 can be defined in 
terms of the inclusive double-differential cc cross section in x and by 

In this paper, the cc cross section is obtained by measuring the D* production cross 
section and employing the hadronisation fraction /(c — ^ D*) to derive the total charm 
cross section. Since only a limited kinematic region is accessible for the measurement of 
D* mesons, a prescription for extrapolating to the full kinematic phase space is needed. 

Since the structure function varies only slowly, it is assumed to be constant within a given 
and y bin. Thus, the measured F^*^ in a bin i is given by 

(r C^) - ^'■""^^'^^^^ ^ ^*^^ F^e (r - O"^) (2) 

0.t,thcoiep U ) 

where cxj are the cross sections in bin i in the measured region of pt{D*) and ri[D*). The 
value of -^2 '^theo calculated from the NLO coefficient functions [S] . The functional form 
of -P'l'^thco used to quote the results for Fg'^'^ at convenient values of Xi and close to 
the centre-of-gravity of the bin. In this calculation, the same parton densities, charm mass 
{rric = 1.35 GeV), and factorisation and renormalisation scales (y^4mfT~Q^) have been 
used as for the HVQDIS calculation of the differential cross sections. The hadronisation 
was performed using the Peterson fragmentation function. 

The beauty contribution was subtracted from the data using the theoretical prediction as 
described in Sectional At low and high x, this fraction is small but it increases with 
increasing and decreasing x. For the lower x point at highest Q^, the contribution 
from beauty production is about 7% of that due to charm production. The contribution 
to the total cross section from Fl'^ calculated using the ZEUS NLO fit is, on average, 
1.3% and at most 4.7% and is taken into account in the extraction of Fg"^. The size of the 
contribution from Fl is similar to that in other PDFs. 
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Cross sections in the measured D* region and in the and y kinematic bins of Table El 
were extrapolated to the full pt{D*) and r]{D*) phase space using HVQDIS. These bins 
correspond to the and x values given in Table |3J where the F2'^ measurements are 
given. Typical extrapolation factors are between 4.7 at low Q^and 1.5 at high Q^, as in 
Table |3] The following uncertainties of the extrapolation were evaluated: 

• using the Aroma fragmentation correction instead of the Peterson fragmentation 
yielded changes of typically less than 10% and not more than 20%. Although these 
values are not very significant compared to the uncertainties in the data, the two 
corrections do produce a noticeable change in the shape of the cross section as a 
function of x. The most significant effects are in the highest x bins for a given Q^; 

• changing the charm mass by ±0.15 GeV consistently in the HVQDIS calculation and 
in the calculation of F2'^ leads to differences in the extrapolation of 5% at low x; the 
value decreases rapidly to higher x; 

• using the upper and lower predictions given by the uncertainty in the ZEUS NLO PDF 
fit, propagated from the experimental uncertainties of the fitted data, to perform the 
extraction of F2'^ gives similar values to the central measurement, with deviations 
typically less than 1%; 

• changing the contribution of beauty events subtracted from the data by ^50''% gave 
an uncertainty of typically 1 — 2% and up to 8% at low x and high Q^. 

These uncertainties were added in quadrature with the experimental systematic uncer- 
tainties when displayed in the figures and are given separately in Table |3] Extrapolating 
the cross sections to the full D* phase space using the CTEQ5F3 proton PDF yielded 
differences compared to the ZEUS NLO QCD fit of less than 5% for Q"^ > 11 GeV^ and 
less than 10% for < 11 GeV\ 

The data are compared in Fig. IHl with the previous measurement P and with the ZEUS 
NLO QCD fit. The two sets of data are consistent^. The prediction describes the data 
well for all and x except for the lowest Q^, where some difference is observed. The 
uncertainty on the theoretical prediction is that from the PDF fit propagated from the 
experimental uncertainties of the fitted data. At the lowest Q^, the uncertainty in the data 
is comparable to the PDF uncertainty shown. This implies that the double-differential 
cross sections given in Table El could be used as an additional constraint on the gluon 
density in the proton. 

^ The first three points of the previous data were measured at = 1.8 GcV^ and not at 2 GcV^, so 
they have been shifted to 2 GeV^ using the ZEUS NLO QCD fit. All other points were measured at 
the same values. 
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The values of Fg^'^ are presented as a function of at fixed values of x and compared 
with the ZEUS NLO QCD fit in Fig. [7| The data rise with increasing Q^, with the rise 
becoming steeper at lower x, demonstrating the property of scaling violation in charm 
production. The data are well described by the prediction. 

Figure IHl shows the ratio F2^/ F2 as a function of x at fixed values of Q^. The values of F2 
used to determine the ratio were taken from the ZEUS NLO QCD fit at the same values 
of and x at which F2'^ is quoted, and are given in Tabled The ratio F2^ jF2 rises from 
10% to 30% as increases and x decreases. 

11 Conclusions 

The production of D* mesons has been measured in deep inelastic scattering at HERA in 
the kinematic region 1.5 < < 1000 GeV^ 0.02 <y< 0.7, 1.5 < Vt{P*) < 15 GeV and 
\'r){D*)\ < 1.5. The data extend the previous analysis to higher and have increased 
precision. 

Predictions from the Aroma MC underestimate, and those from the CASCADE MC over- 
estimate, the measured cross sections. Predictions from NLO QCD are in reasonable 
agreement with the measured cross sections, which show sensitivity to the choice of PDF 
and hence the gluon distribution in the proton. The ZEUS NLO PDF, which was fit to 
recent inclusive DIS data, gives the best description of the D* data. In particular, this is 
seen in the cross-section da / dri{D*) . 

The double-differential cross section in y and has been measured and used to extract 
the open-charm contribution to F2, by using the NLO QCD calculation to extrapolate 
outside the measured pt{D*) and f]{D*) region. Since, at low Q^, the uncertainties of the 
data are comparable to those from the PDF fit, the measured differential cross sections 
in y and should be used in future fits to constrain the gluon density. 
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Table 1: Measured differential cross sections as a function of , x, Pt{D*) and 
ri{D*) for 1.5 < < iqoo GeV^, 0.02 < y < 0.7, 1.5 < Pt{D*) < 15 GeV and 
|?7(D*)| < 1.5. The statistical and systematic uncertainties are shown separately. 
The ratio of the cross sections for e~p and e^p data are also given with statistical 
and systematic uncertainties shown separately. 
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Table 2: Measured differential cross sections as a function of , x, Pt{D*) and 
ri{D*) for 40 < Q'^ < 1000 GeV^ , 0.02 < y < 0.7, l.h < pt{D*) < 15 GeV and 
\ri{D*)\ < 1.5. The statistical and systematic uncertainties are shown separately. 
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Table 3: Measured cross sections in each of the and y bins for 
1.5 < < 1000 GeV^, 0.02 < y < 0.7, 1.5 < priD*) < 15 GeV and 
\r){D*)\ < 1.5. The statistical and systematic uncertainties are shown sepa- 
rately. The prediction for the a^^^^{D*) contribution from HVQDIS, which was 
subtracted from the data in the extraction of is also shown. 
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0.724 


18 


00035 

KJ •kJkJkj KJ'J 


0.293 
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00060 
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500 


0.01200 


0.338 
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+0.029 
-0.092 


+0.021 
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1.57 
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0.03000 


0.180 
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+0.012 
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Table 4: The extracted values of F2'^ at each and x value. The statistical, 
systematic and theoretical uncertainties are shown separately. The values of the 
extrapolation factor used to correct the full pt{D*) and rj{D*) phase space are also 
shown. The value of the proton structure function, F2, from the ZEUS NLO QCD 
fit used to extract the ratio F2'^ / F2, is also given. 
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Figure 1: The distribution of the mass difference, AM = {Mxmrs ~ ^k-k), for 
D* candidates (solid dots). The AM distribution from wrong-charge combinations, 
normalised in the region 0.15 < AM < 0.165 GeV, is shown as the histogram. The 
solid line shows the result of the fit described in the text. The Mk-k distribution for 
the candidates in the range 0.143 < AM < 0.148 GeV is shown as an inset. 
The fit is the sum of a modified Gaussian to describe the signal and a second-order 
polynomial to describe the background. 
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Figure 2: Reconstructed DIS variables for events with D* candidates (after back- 
ground subtraction) for data (points) compared to detector-level Rapgap predic- 
tions (shaded histograms): (a) -(d) show the distributions for 1.5 < < 1000 
GeV^ , while (e)-(h) are the same distributions but for < Q"^ < 1000 GeV^ . All 
histograms are normalised to unit area. 
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Figure 3: Differential D* cross sections, for e~p and e~^p data combined, as a 
function of (a) , (h) x, (c) pt{D*) and (d) rj{D*) compared with MC predictions. 
The inner error bars show the statistical uncertainties and the outer bars show the 
statistical and systematic uncertainties added in quadrature. Predictions from the 
Aroma (dashed line) and Cascade (solid line) MC programs are shown. The 
ratios of the cross sections for e~p and e^p data are also shown beneath each plot. 
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Figure 4: Differential D* cross sections, for e~p and e^p data combined, as a 

function of (a) Q"^ , (h) x, (c) p^iD*) and (d) rj{D*) compared to the NLO QCD 

calculation of HVQDIS. The inner error bars show the statistical uncertainties and 

the outer bars show the statistical and systematic uncertainties added in quadrature. 

Predictions from the ZEUS NLO QCD fit are shown for rric = 1.35 GeV (solid line) 

with its associated uncertainty (shaded band) as discussed in the text. Predictions 

using the CTEQ5F3 PDF (dashed- dotted line) and an alternative hadronisation 

scheme (dotted line) are displayed. The ratios of the cross sections to the central 

HVQDIS prediction are also shown beneath each plot. 
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Figure 5: Differential D* cross sections, for e~p and e~^p data combined, as a 
function of (a) pt{D*) and (b) rj{D*) for > 40 GeV^ . The inner error bars show 
the statistical uncertainties and the outer bars show the statistical and systematic 
uncertainties added in quadrature. Predictions from the ZEUS NLO QCD fit are 
shown for = 1.35 GeV (solid line) with its associated uncertainty (shaded band) 
as discussed in the text. Predictions using the CTEQ5F3 PDF (dashed- dotted line) 
and an alternative hadronisation scheme (dotted line) are displayed. The ratios of 
the cross sections for e~p and e^p data and for e~p and e^p data combined to the 
central HVQDIS prediction are also shown beneath each plot. 
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Figure 6: The measured ¥2^^ at values between 2 and 500 GeV^ as a function of 
X. The current data (solid points) are compared with the previous ZEUS measure- 
ment (open points). The data are shown with statistical uncertainties (inner bars) 
and statistical and systematic uncertainties added in quadrature (outer bars). The 
lower and upper curves show the fit uncertainty propagated from the experimental 
uncertainties of the fitted data. 
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Figure 7: The measured Fg'^ at x values between 0. 00003 and 0. 03 as a function 
of Q^. The data are shown with statistical uncertainties (inner bars) and statistical 
and systematic uncertainties added in quadrature (outer bars). The lower and upper 
curves show the fit uncertainty propagated from the experimental uncertainties of 
the fitted data. 
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Figure 8: The measured ratio F|'^/F2 at values between 2 and 500 GeV^ 
as a function of x. The current data (solid points) are compared with the previous 
ZEUS measurement ( open points). The data are shown with statistical uncertainties 
(inner bars) and statistical and systematic uncertainties added in quadrature (outer 
bars). The lower and upper curves show the fit uncertainty propagated from the 
experimental uncertainties of the fitted data. 
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